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ABSTRACT. We report on new numerical computations of the set of self-contacts in tightly knotted tubes of uniform circular 
cross-section. Such contact sets have been obtained before for the trefoil and figure eight knots by simulated annealing — we 
use constrained gradient-descent to provide new self-contact sets for those and 48 other knot and link types. The minimum 
length of all unit diameter tubes in a given knot or link type is called the ropelengfh of that class of curves. Our computations 
yield improved upper bounds for the ropelengfh of all knots and links with 9 or fewer crossings except the trefoil. 



1. Introduction 

The study of knots as abstract topological objects has inspired a great deal of fascinating mathematics. But knots 
are no less interesting as physical structures tied in flexible ropes and pulled tight. While it is intuitively clear that tight 
knots organize tension and contact forces to bind tightly and resist unravelling, the details of their structure remain 
mysterious. Even today, there is no explicit mathematical description of any tight knot. 

We define the thickness of a space curve to be the supremal (largest) radius of any embedded tubular neighborhood 
of the curve, and the ropelength of a curve to be the quotient of its length and thickness. Ropelength provides a scale- 
invariant way to measure the total flexibility of a given length of rope. It has been established that there is a curve of 
minimum ropelength in each knot and link type L ([6, 5, 3]) and the ropelength of that curve is called the ropelength 
Rop(L) of the knot or link type. These minimum ropelength curves are called tight or ideal knots. In this paper, we 
give some results from our numerical computations of the shapes of tight knots and links. 

Over the past ten years, many authors have found approximate shapes for tight knots and links using numerical 
methods [9, 7, 4, 1]. We follow previous authors in defining a version of ropelength, Rop p , for space polygons and 
optimizing this polygonal ropelength functional over the space of polygons in different knot and link types. But while 
most other others rely on simulated annealing, we introduce the use of constrained gradient descent for ropelength 
optimization. 

This new method has allowed us to significantly expand the scope and accuracy of existing computations of tight 
knots. In particular, it seems that our method has succeeded at the challenging task of resolving the set of self-contacts 
in all knot and link types with nine and fewer crossings (212 types in all). In the process, we have produced a new 
table of upper bounds for the ropelengths of these knot and link types which improves upon previous results. These 
improvements range from 0.05% for the figure eight knot (compared to the bound of [4]) to more than 8.11% for the 
920 knot (compared to the bound of [13]). For links, these seem to be the first upper bounds reported for almost all of 
the link types we consider. 

This dataset is likely to be useful in the study of tight knots, so we provide here an early view of our results. 
This research announcement will be followed by an expanded paper, "Tightening Knots with Constrained Gradient 
Descent", which describes our methods and results in detail. The filesize limitations of the arXiv forced us to truncate 
the data section of this posting, and to use comparatively low-quality image files for the three-dimensional views 
of tight knots and links in Appendix B. A higher-quality view of these images is provided at http : / / www . cs . 
Washington . edu /home s/piatek/contact_t able/. 
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2. Background material 



We start by defining ropelength more precisely. Suppose 7 is a C 1 space curve, parameterized by arclength. The 
maximum diameter of an embedded tube around 7 is controlled by two phenomena: "self-contacts" of the tube formed 
when sections of the curve far away in arclength approach each other in space, and "kinks" formed by points of high 
curvature on 7. These effects are shown below in Figure 1. 
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FIGURE 1 . The thickness of a smooth curve 7 is controlled by curvature (as in the left picture), and 
self-contacts of the tube around 7 (as in the right picture). 




To make this understanding precise, we need a few definitions. We can define the self-distance map of 7 by 

d(M) = ||7(s)-7(*)l|. 

We then have 

Definition 1. The set dcsd(7) of doubly-critical self-distances of 7 is the set of critical points for d(s, t) with s ^ t. 
Taking the partial derivatives of d, we see that (s, t) £ dcsd(7) if and only if 

( 7 (s) - 7 (t), 7 '(*)>= and ( 7 (s) - 7 (t), 7 \t)) = 0. 

Denoting the curvature of 7 at s by k(s), Litherland et al. proved 

Theorem 2.1. [8] The thickness ofj is the minimum of 

■ 1 , . d(s,t) 
mm — — and mm . 

s k(s) (s,t)£dcsd(7) 2 

We can now define the primary object of interest in our computations: 

Definition 2. The self-contact set or strut set of a space curve 7 is the set of (s, t) £ dcsd(7) with || 7 (s) — 7(i)|| = 
2Thi( 7 ). 

The term "strut", borrowed from tensegrity theory, comes from the fact that minimum-length chords in dcsd(7) 
"hold the curve apart from itself" as the knot tightens. 

In our polygonal knot-tightening problem, we will replace the curve 7 with a space polygon V with vertices 
V\,...,vv and edges ei, . . . , ey. To define the polygonal thickness Thi p (V) of V, we will need an idea of curva- 
ture. 

Definition 3. The minimum radius of curvature (or MinRad) of V at Vi is given by the radius of the unique circle 
tangent to both of the edges which meet at Vi and passing through the midpoint of the shorter one. 

If 9i is the turning angle of V at vu then MinRad(wj) is given by 

min{|ej_i|, |e;|} 



MinRad 



2 tan(e*/2) 
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We will also need a definition of doubly-critical self-distances: 
Definition 4. Let dcsd(V) be the set of (p, q) on V with p ^ q which are local minima of the self-distance function. 

We now define a thickness measure for polygons: 
Definition 5. The thickness Thi p (V) of a space polygon V without self-intersections is given by the minimum of 

minMinRadfwi) and min — ^ . 

i (p,g)edcsd(V) 2 

The polygonal ropelength Rop p of V is then the quotient of the length of V and Thi p (V). As expected, when 
polygons V n with increasing numbers of edges are inscribed in a space curve 7, it is known that Rop p (V„) — > Rop(7) 
under some mild geometric hypotheses [10, 12, 13]. We define self-contact sets for polygons like we do for smooth 
curves. 



3. Quantities computed and how the computation was validated 



Our algorithm minimizes the length of V subject to a family of constraints derived from the distance and MinRad 
functions in Definition 5. We report the polygonal ropelengths of the minimized configurations in summary form in 
Appendix A. We also report preliminary computations of the ropelengths of piecewise-smooth curves obtained from 
our polygons by replacing the corners of the polygon with small circular arcs. These provide tentative upper bounds 
on the ropelengths of the knot and link types listed. We also provide self-contact sets for 50 of these knots and links. 



Though our ropelength bounds for knots improve on those of previous authors, this is an uncertain measure of 
their quality. After all, there is no way to check the accuracy of an upper bound for the ropelength of any knot, since 
no exact value for the minimum ropelength of any knot type is known. On the other hand, the minimum ropelength 
is known or explicitly conjectured for some link types [3, 2]. A comparison of our results to some of these known 
examples appears below. 






Link name 
Edges 

Polygon length 
Upper bound 
Smooth length 
Relative error 



Hopf link(2i) 
216 

25.1439 
25.1388 
8?r 

0.02% 



3-link chain (2i#2i) 
384 

41.7131 
41.7086588 

127T + 4 

0.02% 



Borromean rings (6f) 
630 

58.0300 
58.0145 
58.006 
0.01% 



The results above lead us to suspect that most of our smooth ropelength bounds are within 1 to 2 hundreths of a percent 
of the corresponding minimum ropelength values, but we must be cautious. The ropelength "landscape" for knots is 
quite complicated, and we have already discovered a number of local ropelength minima for 8 and 9-crossing knots 
which are very different from the (apparant) global ropelength minima for these knot types. If our gradient descent 
algorithm has been trapped by one of these local minima, the figures we report might be an accurate computation of 
the ropelength of the local minimizer, but far off from the true minimum ropelength value for the knot type. 

We checked our computation of the self-contact sets by comparing our computed self contact set for the Borromean 
rings to the contact set for the ropelength-critical configuration provided by [2] . The results appear in the Figure below. 
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The polygonal configuration of the Borromean rings discovered by our software has total length about 29.01. The 
top plot shows the lower triangle of the square [0, 29.01] x [0, 29.01] representing pairs of arclength values (s,t) on 
these polygonal rings. To describe the position of a point on this 3-component link with a single arclength value, we 
use the convention that arclength values in [0, 9.67) refer to points on the first component, values in [9.67, 19.34) refer 
to points on the second component, and values in [19.34, 29.01) refer to points on the third component of the link. 
These breaks are reinforced by the colored bands running up the diagonal of the plot, which correspond to the colors 
of the different components of our tight links in the 3d renderings of Appendix B. 



The s and t position of points on the plot is indicated by the labels running up the diagonal, where the three special 
values representing breaks between components are lifted away from the other labels. The breaks between components 
are also shown on the plot by the alternating checkerboard pattern of the background colors. 

We then locate every pair (s, t) e dcsd(V) with || V(s) - V(f)|| < Thi p (V) + 10~ 5 . In the context of our gradient 
descent algorithm, these points represent active distance constraints. Our code computes Lagrange multipliers for all 
these active constraints, which represent self-contact forces borne by these tube contacts in the tight knot. We center 
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a dark green box at every such (s, t) value with a nonzero Lagrange multiplier. The size of the box represents the 
average 1 edgelength of the polygon V. We chose this size to represent the expected error in computed self-contact 
positions introduced by approximating a C 1 ropelength-minimizer by the polygon V. 

As we see from the plot, no tube around a component of the link is in contact with itself (so the three cream-colored 
triangles near the diagonal are empty). But each of the components makes contact with the other two, as shown by 
the boxes plotted in the purple and cream-colored rectangles forming the remainder of the plot. We can see that the 
contacts break up naturally into "lantern-shaped" structures. 

This link has been studied by Cantarella, Fu, Kusner, Sullivan, and Wrinkle, who provide a ropelength-critical 
configuration in [2]. In the bottom plot, we compare one "lantern" formed by 608 of these boxes to the self-contact 
set predicted by these authors, which is represented by a red line. In this plot, the arclength distances labelled on the 
s and t axes do not correspond to a region of the plot above, but merely indicate the scale of the plot. We can see that 
the agreement between theory and computation is generally within one edgelength. 

Appendix B contains similar plots of our computed self-contact sets for 50 knots and links from our collection of 
minimized examples. To the left of each self-contact plot, we provide a 3d rendering of the corresponding tight shape. 
For links, the diagonal and right-hand side of the triangular self-contact plot are colored gray, orange and green in 
correspondance with the colors of the components in the rendering. The start of each component in the rendering 
is denoted by the tube coming to a point. The black bands on the tubes correspond to the arclength tick marks on 
the plot at right. The table also contains the polygonal ropelength of the knot (top number, slightly higher) and the 
corresponding smooth ropelength upper bound (bottom number, slightly lower), as well as the number of edges in the 
configuration plotted. 



4. Conclusions 

The major contribution of Appendix A is the provision of ropelength figures for links. This allows us to check the 
accuracy of numerical ropelength minimizations against theoretical results for the first time. We are happy to report 
that our method passes this test for the cases we examined. 

The pictures in Appendix B are considerably more evocative. It is evident from first inspection that the contact sets 
of tight knots and links seem to contain a fairly small number of commonly repeated patterns. Some of these, such as 
the "steps" pattern first seen in the trefoil knot, change shape from knot to knot. But others (such as the "winged" pat- 
tern in the figure eight knot or the "lantern" shape seem in the Borromean rings) seem to remain remarkably consistent 
throughout our computations. The reader may notice many other examples as well. Isolating and understanding some 
of these structures could provide us with a "construction kit" for tight knots, much like the analysis of the simple chain 
in [3] led to the construction of an infinite family of tight links. 
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Our final polygons are almost-equilateral, so this is a good approximation of the lengths of the edges incident to the pair (s, t). 
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Appendix A. Table of polygonal ropelengths and ropelength upper bounds 



The table shows new polygonal ropelengths and ropelength upper bounds for 212 knots and links with 9 or fewer 
crossings. From left to right, there are four columns: the name of the knot or link, the polygonal ropelength Rop p , 
the corresponding ropelength upper bound Rop, and the previous best ropelength upper bound we could find in the 
literature together with the percentage improvement in ropelength. For these last figures, we used the papers [1,4, 11]. 



Link 


Rop p 


Rop 


Previous 


2\ 


25.1439 


25.1388 




3i 


32.7490 


32.7448 


32.7433864 (-%) 


4i 


42.0997 


42.0928 


42.1158845 (0.05%) 


4? 


40.0247 


40.0169 




5i 


47.2156 


47.2016 


47.51 (0.64%) 


59 


49.4840 


49.4704 


49 73 f0 52%) 


r- 9 

5? 


49.7874 


49.7723 




Dl 


^fi 7^1 fi 


^fi 71 ^.n 


o / .11 i^u.uy /o j 


R 
t>2 


C7 A/1 ^1 


C7 0971 


c7 /M 71 Q/\ 


t>3 


^7 8fi n 9 




Oo.4o ^l.Uo/oJ 


6? 


54.4068 


54.3893 




61 


56.7132 


56.7028 




b 3 


1 1 A1 
Do. 11D1 


Do. 1U44 




6? 


57.8334 


57.8170 




61 

Z 


58.0300 


58.0145 




o 


50.5865 


50.5745 




7i 


61.4319 


61.4109 


61.89 (0.77%) 


7 2 


63.9165 


63.8956 


65.36 (2.24%) 


73 


63.9539 


63.9327 


64.35 (0.64%) 


7 4 


64.2960 


64.2724 


65.63 (2.06%) 


7 5 


65.2802 


65.2609 


65.70 (0.66%) 


7 ' r- 

< 6 


fi^l 71 


fi^i 701 9 


fifi 1 7 fn 7%^ 


' 7 


65 6316 


65 6108 

U V . U J- U U 


66 09 (0 72%") 


7? 


64.2585 


64.2353 






65.0467 


65.0274 




n 


65.3743 


65.3561 




7 2 

' 4 


65.0971 


65.0759 




'5 


66.2400 


66.2186 




7 2 
'6 


66.3494 


66.3372 




7 2 
' 7 


55.5451 


55.5311 




7| 


57.8043 


57.7948 




7? 


65.8275 


65.8090 




8i 


71.0484 


71.0241 


71.44 (0.58%) 


8 2 


71.4327 


71.4107 


71.91 (0.69%) 



Link 


Ro Pp 


Rop 


Previous 




83 


71 


1880 


71 


1655 


71 


56 


(0 


55%) 


84 


72 


0301 


72 


0049 


72 


41 


(0 


55%) 


85 


72 


2100 


72 


1878 


72 


70 


(0 


7%) 


8 6 


72 


5005 


72 


4791 


72 


93 


(0 


61%) 


87 


( L 


0/1/17 
Z44 t 


70 


OOfl/1 
ZzU4 


70 
/ Z 


Do 


(C\ 
[\J 


56%) 


8s 


73 


3533 


73 


3334 


73 


88 


(0 


73%) 


89 


72 


4717 


72 


4461 


72 


96 


(0 


7%) 


810 


73 


4279 


73 


4095 


73 


86 


(0 


6%) 


811 


73 


5029 


73 


4802 


76 


70 


(4 


19%) 


812 


74 


0291 


74 


0098 


74 


61 


(0 


8%) 


813 


72 


8291 


72 


8045 


73 


29 


(0 


66%) 


814 


73 


9226 


73 


8991 


74 


93 


(1 


37%) 


815 


74 


3344 


74 


3134 


74 


82 


(0 


67%) 


816 


74 


9213 


74 


8962 


75 


47 


(0 


76%) 


817 


74 


5276 


74 


5071 


75 


08 


(0 


76%) 


818 


74 


9420 


74 


9252 


75 


44 


(0 


68%) 


819 


61 


0734 


61 


0430 


61 


35 


(0 


5%) 


820 


63 


1530 


63 


1146 


64 


11 


(1 


55%) 


821 


65 


5504 


65 


5298 


65 


91 


(0 


57%) 


8? 


68 


5198 


68 


4884 










^9 

82 


1 1 




1 1 


loo 1 










8? 

°3 


72 


7498 


72 


7291 










81 


72 
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72 


5908 










8g 


74 


0039 


73 


9826 










81 


73 


2932 


73 


2502 










8? 


74 


4165 


74 


3885 










81 


73 


7849 


73 


7702 










81 


74 


0620 


74 


0386 










8 io 


73 


6890 


73 


6684 










8?i 


73 


0115 


72 


9899 










812 


74 


0194 


73 


9140 










8 13 


74 


1685 


74 


1501 










8 14 


73 


7000 


73 


6775 










8? 5 


64 


3305 


64 


3086 










8l6 


66 


8434 


66 


8315 










8? 


72 


2883 


72 


2649 










83 


72 


9544 


72 


9360 










8§ 


74 


9366 


74 


9139 










81 


77 


8544 


77 


8314 










81 


73 


4286 


73 


4061 










8| 


74 


7680 


74 


7468 










8? 


60 


6065 


60 


5888 
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Link Rop p Rop Previous 



Link Rop p Rop Previous 



8| 65.0637 65, 

8^ 70.1904 70. 

8 3 w 68.9823 68. 

8f 75.2901 75, 

8| 67.4772 67. 

8^ 66.4140 66. 

91 75.7507 75, 

9 2 79.3059 79. 

9 3 78.5867 78, 

9 4 78.4117 78, 

9 5 79.7877 79, 

9 6 80.1326 80, 

9 7 80.6597 80, 

9 8 80.5651 80, 

9 9 79.9574 79, 
9i 79.8161 79, 
9n 80.3650 80, 

912 80.1275 80, 

91 3 80.6525 80, 

914 80.1512 80, 

91 5 82.1665 82, 

91 6 80.1446 80, 

91 7 80.5792 80, 

91 8 81.6230 81, 

91 9 82.1828 82, 

9 20 80.2543 80, 

9 2 1 81.1336 81, 

9 22 81.0714 81, 

9 23 81.3164 81, 

9 24 80.99 3 3 80, 

9 25 81.1873 81, 

9 26 80.9352 80, 

9 27 81.9473 81, 

9 28 81.5694 81, 

9 29 81.8708 81, 

9 30 81.8 5 6 7 81, 

9 3 1 81.7436 81, 

9 32 81.5775 81, 

9 33 82.8451 82, 

9 34 82.32 1 3 82, 

9 35 79.2495 79, 

9 36 81.0579 81, 

9 37 81.5845 81, 

9 38 81.8119 81, 

9 39 81.9490 81, 

9 40 81.7008 81, 
9 44 81.4929 81, 

9 42 69.6 1 33 69, 

9 43 71.70 6 2 71, 
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1460 


r\9 

9| 


81 


0999 


81 


0723 


r\9 

9g 


80 


3391 


80 


3181 


r\9 

9?o 


80 


3964 


80 


3693 


n9 

9?i 


82 


0513 


82 


0271 


y i2 


81 


9983 


81 


9738 


r»9 

9? 3 


79 


3586 


79 


3319 


r»9 

9? 4 


80 


7420 


80 


7171 


9? 5 


80 


5805 


80 


5550 


9? 6 


81 


4190 


81 


3966 


rt9 

9? 7 


81 


2589 


81 


2307 


9? 8 


82 


2739 


82 


2429 


y i9 


79 


4899 


79 


4628 


o2 
^20 


80 


4176 


80 


3956 


9ii 


81 


2545 


81 


2356 


9i 2 


81 


1318 


81 


1058 


9i 3 


80 


4552 


80 


4290 


9i 4 


82 


5780 


82 


5207 


9i 5 


81 


8205 


81 


7924 


9le 


82 


1200 


82 


1015 


r>9 

9i 7 


81 


3442 


81 


2046 


r»9 

9l 8 


81 


3869 


81 


3740 


r»9 


82 


1780 


82 


1567 


9| 


82 


2789 


82 


2517 


9§i 


80 


6125 


80 


5936 


/i9 

9 32 


81 


4352 


81 


4122 


r\9 

9 33 


82 


2167 


82 


1951 


r»9 

9i 4 


81 


8743 


81 


8552 


r\9 

9§ 5 


81 


3504 


81 


3284 


9§6 


80 


7280 


80 


7103 


9 37 


81 


9353 


D 1 

81 


9101 


q2 


82 


7480 


82 


7206 


q2 

y 39 


81 


9251 


81 


8973 


q2 

y 40 


82 


0090 


81 


9822 


q2 

y 41 


83 


6167 


83 


5979 


q2 

y 42 


83 


6493 


83 


6292 


q2 

y 43 


66 


3264 


66 


3049 


q2 

y 44 


72 


2518 


72 


2255 



72.23 (0.82%) 
75.51 (0.81%) 
69.35 (1.02%) 
75.61 (0.69%) 
74.94 (0.92%) 
74.50 (0.65%) 
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Link Rop p Rop Previous 



91 5 71.3593 71.3490 

91 6 73.9646 73.9377 
%\ 7 69.9602 69.9432 
9i 8 73.6781 73.6545 
9% 66.0806 66.0658 
9§ 69.3690 69.3483 
9§! 70.5942 70.5699 
9| 2 72.9916 72.9685 
9| 3 68.0369 68.0305 

9 2 54 71.0385 71.0185 

9 2 55 73.8423 73.8217 
9| 6 75.2439 75.2245 
9| 7 73.8408 73.8194 
9§ 8 74.1911 74.1697 
9^ 9 73.0481 73.0305 
9§ 73.5734 73.5553 
9§! 69.3978 69.3840 

9? 81.2897 81.2323 

9\ 82.4507 82.4004 

9 3 3 82.3127 82.2861 

9| 82.5449 82.5208 

9^ 80.7974 80.7714 

9% 81.0235 81.0054 

9 7 82.1986 82.1535 

91 81.1631 81.1408 

9| 81.6200 81.5735 

9? 82.3446 82.3259 

9^ 82.0323 82.0137 

9 3 12 82.6345 82.5740 

9? 3 72.2098 72.2008 

9? 4 74.5697 74.5492 

9? 5 74.3877 74.3655 

9? 6 75.0664 75.0430 

9 3 17 74.2972 74.2779 

9? 8 72.5059 72.4741 

9? 9 72.7143 72.6859 

9^ 76.3557 76.1829 

9| x 74.9369 74.9212 

9f 85.5620 85.5115 
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Appendix B. Table of Self-Contact Sets 
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